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c<  absorptivity 
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INTRODUCTION 

Today,  more  than  ever  before,  the  role  of  accurate  error  analysis 
in  approxinate  procedures  has  become  e:ctremely  vital.  With  the  advent 
of  high  speed  digit^d  conr^uters,  the  solution  of  ir.any  difficult  and  soins- 
tiroes  complex  problems  has  been  made  possible  using  approximation  tech- 
niques. The  number  of  appi^ximation  techniques  available  are  many.  Un- 
fortunately one  does  not  have  a  definite  idea  of  the  error  which  is 
present  using  such  techniques,   -■/  „.  ^-  •    v   • 

jji  this  report,  concerned  '.■Jith  radiation  heat  transfer,  the  purpose 
is  to  formulate  a  procedure  to  determine  bounds  for  the  radiant  flux  along 
the  surface  between  tvjo  parallel  plates  of  finite  extent.  The  results  are 
compared  vdth  exact  and  other  approximate  methods.  Once  good  upper  and 
lower  bounds  have  been  obtained,  the  maximum  possible  error  can  be  deter- 
mined. Emphasis  should  be  placed  on  the  ability  of  this  method  to  yield 
an  approximate  solution  of  known  accuracy. 

The  bounding  principle  is  a  method  in  which  converg^it  upper  and 
lower  bounds  are  obtained  in  order  tc  approximate  the  solution  to  the 
problem.  The  formulation  of  the  upper  and  lower  bounds  is  established  by 
bounding  functions.  These  bounding  functions  essentially  are  intermediate 
problems  to  the  actual  problem.  By  varying  these  functions,  the  solution 
of  the  intermediate  problem  is  always  numerically  greater  than  the  exact 
solution  in  one  case  and  always  numerically  lower  than  the  exact  solution 
in  the  other  case.  Applying  this  procedure  systematically,  the  exact 
solution  is  then  bounded  above  by  a  solution  which  is  greater  and  bounded 
below  by  a  solution  which  is  less. 


The  method  is  practical  in  application  since  it  is  not  necessary  to 
obtain  complete  convergenceo  Once  good  bounds  are  established,  it  is  not 
onlj'-  possible  to  obtain  a  good  approximate  solution,  but  also  to  specify 
\d.th  certainty  the  raaximum  possible  error.  Thus,  the  difficult  problem 
of  error  analysis  is  avoided. 

The  procedure  for  vjhich  the  bounding  principle  is  established  is  one 
which  has  been  developed  for  determining  bounds  for  the  solution  of  differ- 
ential equations.   H  5  I^l  *  [3J 

A  radiation  problem  v;as  chosen  for  th^  pirrpose  of  demonstrating  the 
application  of  the  bounding  principles  to  integral  equations  which  arise 
in  radiation  heat  transfer  problems.  Since  the  solution  of  integral 
equations  in  general  is  a  difficult  task  even  for  those  who  are  familiar 
with  them,  a  simple  physical  situation  was  selected  for  this  problem.  In 
this  way,  the  essential  features  of  the-  formulation  are  retained  without 
the  difficulties  which  woxild  be  involved  in  a  complex  geometrical  arrange- 


ment 


LJ      Kuffjbers  in  brackets  designate  References  at  end  of  report. 


DESCRIPTICN  OF  RADIMT  FLUX  PROBLEM 

The  problem  is  to  detemiine  the  rate  of  the  combined  radiation  flux 
along  the  surface  between  two  finite,  non-black,  parallel  plates. 

Consider  the  figure  below  in  which  a  sketch  of  the  parallel  plates 
is  sho\TO, 


k\^^\^^^^^^ 


h 


^^^^^=^^^^^^ 


•r- 


Fig.  !♦  Two  Non-Black  Parallel  Platec. 

Each  plate  is  of  length  L.  The  plates  are  directly  opposed  and 
extend  indefinitely  in  the  direction  nornial  to  the  plane  of  the  figure. 
A  distance  h  separates  the  two  parallel  plates.  Each  plate  is  considered 
as  a  grgy-body  radiator  \chich  both  radiates  and  reflects  in  a  diffuse 
nanner.  The  gray-bod?/-  eaissivities  of  both  plates  are  the  same.  Thus, 


6  =  oc 


/-p 


(1) 


The  ffisdium  separating  the  two  plates  may  contain  either  a  non- 
participating  gas  or  a  vacuum. 

The  surface  tesraera-iure  of  the  plate  is  a  function  of  the  distance 
c2.ong  the  plate.  For  this  problem  tho  teimserature  is  assumed  constant 
and  both  plates  aro  ntrhitainad  at  tho  aame  uniforra  tetrperature. 

The  net  heat  transfer  from  the  system  vdJLL  be  due  to  the  radiant 
energy  which  is  lost  through  the  ends  of  the  gap  to  the  environment, 

Tne  rate  of  the  combined  radiation  flux  (omitted  and  reflected) 
leaving  a  unit  area  dA.,  at  x  is  designated  as  B(x)«  The  emitted  part 
of  the  flux  is  given  by  the  well-knovni  relation,  e<rT^,  The  coiount  of 
the  reflected  portion  of  the  flxix  is  PE(x),  where  H(x)  represents  the 
radiant  energy  arriving  at  x  per  unit  time  and  area.  Thus  the  rate  of 
the  coTiibined  radiation  flux  can  be  expressed  as 

B(.X)    =  ^'(rJs^'^)    ■*■  ^HC7^)  (2) 

(The  B,  H  notation  follows  the  notation^ used  in  illumination  engineering 
and  is  the  notation  used  by  Sckert,  [UJ ) 

Writing  an  energy  balance  for  an  elemental  area  of  surface  leads  to 

the  integral  equation  for  the  radiant  fliix  between  two  non-black,  parallel 

plates,  [Uj  ,  [5]  ,  [6]  The  integral  equation  is  as  follosfs 


(3) 


It  may  be  noted  here  that  the  method  used  in  most  standard  textbooks 
to  determine  the  reflected  portion  of  the  radiant  flux  makes  the  assump- 


tior.  that  it  is  constant  for  all  values  of  x  rather  than  being  a  function 
of  X. 

The  integral  equation  is  transformed  into  a  dinxensionless  form  by 
dividing  through  byeff-T^  ,  whore  Tq  is  sovjs   arbitrarily  chosen  reference 
tezperaturo.  Thus,  equation  (3)  becomes 


.-^  .  ^  r-2  ^   ,  — ^ 


Define  the  following  expressions: 


ecr  L 


^(Y)  ^  _Sjui_  (5) 


e^ 


(6) 


Js(X)    =   Ts  (z) 

z 

X  and  Y  are  diniensionless  coordinates  defined  as  X  =  x  and  1  "»  y  , 
respectively.  Substituting  these  quantities  into  equation  (3a)  gives 


^U)    =  JsU)     .  ^f    ^^1^ cl^-L      (3,) 


~i 


[iy-t-^ir-''T 


The  denominator  of  the  integral  expression  becomes,  after  factoring  an 
L/2  out  of  the  expression 

which  yields 


my-4H:^] 


\V2.f 


Define  H  as  the  gap  spacing  ratio  H  -  h   and  substitute  into  equation 
(3b).  This  yields 


METHOD  OF  DETERMINING  BOUNDS 

Vpper   and  lower  bounds  are  established  in  terms  of  bovmding  functions. 
These  bounding  functions  are  essentially  intermediate  problems  to  the  actxial 
problem.  The  solution  of  the  intermediate  problem  is  always  numerically 
greater  than  the  exact  solution  in  one  case  and  always  numsricaLly  lower 
than  the  exact  solution  in  the  other  case.  Systematic  procedures  can  then 
be  developed  to  provide  for  convergence  of  these  bounding  functions  to  the 
exact  solution,  ,  , 

Let /^  (X)  represent  a  continuous  function  in  the  interval  0  =•  X  »  1, 
If  ^  (X)  is  substituted  for  fl(X)  in  the  integral  equation  (3c),  the  equation 
will  not  be  satisfied  in  general.  Therefore,  ^■^^(X)  is  not,  in  general,  the 
solution  of  the  combined  radiation  flux  problem  with  surface  temperature, 
j7g(X),  However,  /5*(X)  does  represent  the  exact  solution  of  a  problem  corre- 
sponding to  some  other  temperature.  This  will  be  denoted  byJ7*(X),  Using 
equation  (3c),:7*(X)  can  be  determined  explicitly  by  direct  substitution  of 
^*(X)  for|5(X),  Thus 


J./  ,' 

where  U  is  defined  as  U  -  (X  -  X)^  ♦  H^.  it  becomes  apparent  from  equation 
(7)  that  for  a  given  solution,  the  determination  of  the  exact  problem  is  a 
straight  forward  procedure. 
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Suppose  that  a  particular  jS*(X)  «  i6  (X)  can  be  chosen  such  that  the 
corresponding  tenperature  j'^Oi)   ■  ^(X)  satisfies  the  following  inequality: 

J;,(X)   ^    JuiX)  O^X^  /  (8) 

The  net  radiant  flux  distribution  3  (X)   corresponds  to  a  situation 
where  the  surface  temperature  is  greater  than  the  desired  temperature, 
C7g(X).  Intuitively,  one  would  conclude  that  the  effect  of  an  increase  in 
the  temperature  would  caiise  a  corresponding  increase  in  the  net  radiant  flux* 
This  can  be  expressed  as 


/^(X)  =   s.(x) 


O^X^I  (9) 


Similarily  if  ^*(X)  -  ^(X)   can  be  chosen  such  that  J*(X)  -  J  (X) 


satisfies 


J^  (X)   ^    Js  (X)  O^X^  I        (10) 


then  it  foUotfs  that 


6,  (X)  ^  ^(X)  OM  X5  /    (11) 


From  equations  (9)  and  (11),  upper  and  lower  bounds  for  &{X)   are  th\xs 
established. 


^j.  (X)  ^   ^  (X)  ^  I5u  (X)  O^  X^  I      (12) 


The  problem  now  becomes  one  of  determining  a  practical  procedure  to 
obtain  good  upper  and  lower  bounds  in  such  a  manner  that  the  deviation 
is  sufficiently  small.  The  answer  is  to  begin  by  choosing  a  /5*(X)  =  y^(X) 
such  that  two  requirements  are  satisfied.  One,  JA^)   should  be  approximately- 
equal  to  J  (X)  and,  two,  J^CX)  -  Cr^(X).  Likewise,  the  choice  ^*(X)  =  /^(X) 
must  be  made  such  that  CZ  (X)  as^  (X)  subject  to  the  condition,  C^  (X)  =  ^(X). 

Is  X        i> 

Moreover,  emphasis  should  be  placed  upon  the  fact  that  to  obtain  convergent 
bounds  not  only  laust  J7  (X)  and  C7  (x)  approach  J'(X),  but  the  inequalities  of 

XL  X  S 

equations  (8)  and  (10)  must  be  satisfied  during  the  limiting  process.  It  is 
apparent  that  the  choice  of  ^*(X)  must  result  in  a  one-sided  convergence  for 
each  bound. 

The  first  step  in  determining  bounds  is  to  assume  a  series  form  of  the 
solution  for  3*(X),  Due  to  the  symmetry  of  the  problem,  assume  a  power 
series  solution  for  /S*(X)  and  fi*(I)  in  equation  (7)  to  be 


^''(X)    =Y.    ^'^^'^  (13) 


and 


M 


Substitution  into  equation  (7)  leads  to 


^"(k;  -^  A^Y""  (lu) 


j^         w  M 


y-/ 


J'"(X)      =  £  ^y"  -^JLY^   ^  [ll  dY  (7a) 


10 

The  definite  integrals  are  now  defined  as 

z« 

J-,  u"\ 

(15) 

i 

and  the  indefinite  integrals  . 

are  written  as 

i 

•J 

1  ^% 

(16)        I 

Prota  a  table  of  integrals  \j\ 

,  an  equation  is  written  for  each  value 

of             i 

■! 

N  in  equation  (l6)  for  N  «  0, 

1,  2.  A  general  equation  is  then 

developed 

for  any  value  greater  than  2, 

The  evaluation  of  the  indefinite 

integrals 

1 

is  given  in  Appendix  A,  The 

equations  for  Cjj  are  as  follows; 

i 

LaJO     ~          T-~X. 

i 
(l6a) 

^    .■-     '           i    '    '' 

'1 
(16b) 

■•  -  Sj 

•■'■•■   -.'^ 

i 
i 

u 


C^^   =  /^(/-X^l/'")    ^  /'9'-X^-X^'-  ^y  (16c) 


t- 


iA/-Z)U'^^  A/-Z  ~  A/ -2. 


Evaluating  the  definite  integral  in  equation  (l5)  for  N  -  0  at  its  limits 
of  -1  and  ♦l  gives 


.iJO 


=   C^o(xj.o)  -  C^oU.'i.o) 


^^o    =  '-X         ^       /^/  (17) 

where  U^^  -  (1  -  X)^  +  H^  and  U_^  -  (l  +  X)2  +  H^.  For  N  »  1,  the  equation 
is 


2^,  =  C^/^./-^/-  C^,(x.-io) 


Z^,  =   /-/-//^  ^  Z-^/^^//'  (18) 


For  N  "  2 


^^,  =  (f-^z  ("/,  10 )  -  <^v,  r/.  -i.o) 
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^>^/-^  =   /n(i-X^ul^ )   -    Jn   (f-X^u!^) 


^  /'-X'-XH^-H'-    -t  X!±x!±XJdl:j£_      (19) 


■^1  -' 


^/a. 


Finally,  N  >  2  leads  to 


^^>z  =    C^>:.  (X,  /.o)  -  C^>^  CX,-l.O) 


/       -  i-l)  (20) 

Exaniining  equation  (7a),  the  next  step  becomes  one  of  determining  the 
coefficients,  A^  ,  from  the  assiuned  power  series  solution  for  ^"^(X),  This 
is  accon5>lished  by  using  the  method  of  collocation. 

The  method  of  collocation  is  a  relatively  simple  and  straight  forward 

approximation  method  which  in  essence  is  a  weighted  residual  procedure.  It 

is  used  here  to  "force"  J^(X)   and  J,  (X)  to  be  approximately  equal  to  J^(X) 

s 

in  the  interval  0  »  X  »  1,  It  is  not  satisfactory,  in  general,  to  let  J'(X) 
and  ^(X)  approach  C7(X)  in  the  usual  manner.  In  order  that  upper  and  lower 
bounds  may  be  established,  the  requirement  that  the  inequalities  of  equations 
(8)  and  (10)  must  also  be  maintained.  It  then  becomes  convenient,  therefore. 
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to  let  J  (X)  and  jAl)   approach  a  new  value,  J7(X)  +6,  where  5  is  a 

parameter.  By  suitable  choices  of  S,  the  inequalities  can  be  maintained. 

Also,  it  is  convenient  to  define  a  residual  function,  R  _(X),  as 

No 

/C,r  (^^   =  J^'^C^)  -  Jf(^)  (21) 


The  method  of  collocation  consists  of  selecting  a  number  of  equally  spaced 
points  for  X  in  the  interval  0  =  X  «  1  and  equating  the  values  of  the 
residuals  to  S  at  each  one  of  these  points.  This  leads  to  a  system  of 
simultaneous  equations. 


/?^s(^i^     -  /  =   ^ 


A/^  {^/^  ~  ^   -   ^  (22) 

where  i  »  1,  2,  .  ,  .  ,  (M  ♦  1),  The  A  «s  can  then  be  determined  from 

2  ^ 

any  matrix  solution  program. 

Once  the  coefficients  have  been  obtained,  the  next  step  in  the  bounding 
technique  is  to  consider  the  residual  equation  (21)  again.  To  insure  bounds, 
the  residuals  must  be  checked  to  determine  if  the  proper  inequality  has  been 
satisfied  at  all  points  in  the  interval  0  =  X  =  1,  For  each  value  of  6, 
the  residuals  are  examined  at  101  points  to  determine  the  approiximate  maximum 
and  minimum  values  of  the  residual. 

An  upper  bound  for  the  net  radiant  flux  is  obtained  by  choosing  a  suf- 
ficiently large  positive  value  of  S,  The  residuals  and  the  A's  are  then 
calculated.  If  the  minimum  residual  is  positive,  then  the  inequality  (8) 
is  satisfied  and  the  A's  can  be  substituted  into  equation  (?)•  The  resulting 


Ii4 


6  (X)  is  an  upper  bound  for  ^iX)»     If  the  bound  scheme  is  repeated  for 
successively  smaller  values  of  S  and  the  rainimuin  residual  approaches  zero 
but  still  remains  positive,  a  least  upper  boxmd  will  be  attained  for  |5(2)» 

The  lower  bound  for  the  radiant  flux  is  similarly  obtained  by  choosing 
a  sufficiently  large  negative  value  of  8  and  determining  the  residuals  and 
the  A's.  If  the  maximum  residual  is  negative,  then  the  inequality  (10)  is 
maintained  and  the  A's  can  be  substituted  into  equation  (7).  The  resulting 
/^l(X)  is  a  lower  bo\ind  for  y^(X),  Better  lower  bounds  can  be  obtained  by 
letting  6  approach  zero  and  verifying  that  the  sign  of  the  raaxiraTm  residual 
remains  negative, 

A  schematic  for  illustrating  how  the  maximum  and  minimum  residuals 
vary  with  different  values  of  6  is  shown. 


Greatest  Lower  Bound 


Maximum  Residual 


Minimum  Residual 


•Least  Upper  Bound 


Fig.  2,  Maximum  and  Minimum  Residuals  versus  S 


f 


^■ 


Based  upon  the  definition  of  the  residual,  it  can  be  seen  that  for  the  least 
upper  bound  and  the  greatest  lower  bound,  the  inequalities  of  equations  (3) 
and  (10)  are  maintained. 

The  best  approximation  to  the  exact  solution  is  the  mean  value  of  the 
upper  and  lower  boxindlng  solutions.  Thus 

L  (/)   =   A  (/)    -^    /z  U)  (23) 

is  the  mean  value.  The  maximum  percent  error  associated  id.th  this  mean 
value  is  expressed  as 


Maximum  Percent  Error       =  /^^(^)~  $l(X)     (iqq\  ^^li) 

f>AX} 


NUMERICAL  RESULTS 

The  numerical  results  obtained  for  the  net  radiant  flux  along  the  surface 
of  the  plate  are  compared  vrLth  a  solution  obtained  by  direct  iteration  on  an 
electronic  computer.  The  iterative  solution  is  given  by  Sparrow.  (jSJ  This 
solution  is  used  as  a  standard.  In  all  of  the  cases  considered,  "OjX)   **  1»0 
to  agree  with  the  work  in  [6]  ,  The  numerical  results,  corresponding  to  three 
cases,  are  given  in  the  Tables  I,  II  and  III.  Table  IV  shows  the  coefficients 
of  the  bounding  functions  for  the  three  cases  considered.  The  S  ' s  selected  in 
the  three  cases  are  the  best  values  for  obtaining  good  upper  and  lower  bounds. 
Figure  3  shows  a  comparison  of  the  mean  value  and  the  exact  solution  for  a 
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particular  case.  Figure  U  shows  how  the  bounds  vary  along  the  surface  of  the 
plate  for  different  values  of  S« 

The  numerical  work  was  done  on  an  IBM  ll;10  digital  computer  using  8 
place  arithmetic. 

Five  terms  were  selected  in  the  series  solution  for&*(X)  and  ^{1) 
in  equations  (13)  and  (lii). 
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Table  17 
Coefficients  of  the  Bounding  Functions 

Collocation  N  =  5j  H  «=  2.0;  (>  =  0,9;  D^(X)   =  1.0 


Lower  Bound  ,        Upper  Bound 
(  6  =  -0.1  X  10"^)    (  S  =  0.1  X  10--*) 


h 

0.16696929  X  lO-'- 

0,16697050  X  10- 

^ 

-0.159U0978 

-0.l59i;Oiai 

0.2l89i;532  x  10'^ 

0.21851796  X  10-1 

-0.97139^90  X  10-3 

-0.887929li3  x  10"^ 

^ 

-0.21857891  X  10"^ 

-0.26551958  X  10-3 

Collocation  N  =  5;  H  =  1.0;   p  -  0.9;  ZJ  (X)  »  1.0 


Lov;er  Bound         Uoper  Bound 
(  S   -  -0.3  X  10-3)    (  5  «*  0.3  X  10-3) 


^ 

0.25072739  X  10^ 

0.25077759  X  10^ 

^ 

-O.6U818792 

-0.61;83178l 

s 

0.1536371;5  x  lO'^ 

0.15368850  X  10"! 

\ 

0.61335306  X  10-1 

O.613U3168  X  10-1 

^ 

-0.11606111  X  10"! 

-0,11606003  X  10"! 
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Table  IV     (continiied) 


Collocation    N  -  ^j     H  »  0.2j      P  =  0.9j     C7.(X)  ==  1,0 


Lower  Bound    ^  Upper  Bound 

(  5  -  -0.725  X  10-^)       (  5  =  0.2  X  10-^) 

A  O.70U09250  X  10-^  O.7606U76I  X  10^ 

Ag  -0.17770760  X  10^  -0.19198206  X  10^ 

A  -0.38609392  -0.1il710001 

^  2  -1 

A  -O.9U682988  X  lO-'^  -0.102la78U  X  10  -^ 

1; 

A  -0.1i;023075  X  10^  -0.l5lii936l  x  10^ 
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1.70 


1.63 


1*66 


Msaa  Value,  y^aCX) 


Exact* 


1.61itl^  "  "- 


•i 

o 
o 


1.62 


lo60 


1.58  — 


1.56 


1.51; 


1.52 


1.50 


Collocaticxi 


N=5 
H  »  2,0 
^  »  0.9 

J^gCX)    »   1.0 

3"t=  -0,1  X  10 


<r^i"  0.1  X  10 


^  Exact  solution  is  based  on  a  direct  iteration. M  ^ 

I      i      I      I      I 


I      I      ■      ■      t 

0.6       0»8       IcO 


0.2        0.U 

DiDisnsionless  Position  X 
Fig,  3,     Coxparison  of  Mean  Valua  and  Exact*  Solution 
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0.2  OoU  0,6 

Dimensionless  Position  X 
Fig.  k»       CoHparison  cf  Botinds  versvis  S» 
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DISCUSSION  OF  RESULTS 

Th3  present  results  indicate  that  the  convergence  of  the  upper  and  loiter 
bounds  is  obtained  for  the  combined  radiation  f lu:c  along  the  surface  of  the 
plate.  Bounds  are  established  for  each  of  the  three  cases  considered.  The 
percent  error  remains  relatively  constant  along  the  surface  of  the  plate  for 
each  case.  This  indicates  that  the  deviation  between  the  bounds  is  constant 
along  the  surface. 

As  the  gap  sp^-^ing  ratio,  H,  is  decreased  from  2,0  to  0,2,  the  ma::imura 
percent  error  becociss  larger.  Thus,  it  is  concluded  that  better  bounds  are 
obtained  for  larger  gap  spacing  ratios, 

A  comparison  of  the  results  obtained  and  the  exact  solution  shovjs  that 
they  are  not  in  agreemsnt.  It  can  be  seen  in  Fig,  3  that  the  shape  of  the 
curve  for  both  solutions  is  approximately  the  sasne.  The  solu\iions  appear  to 
diffsr  by  a  constant.  This  implies  that  the  disagreement  is  not  too  serious 
and  is  possibly  due  to  either  a  systematic  or  numerical  error.  Further  in- 
vestiration  is  recommended  to  determine  the  nature  of  the  disagreemant. 

In  the  first  case  considered,  it  is  shown  in  Fig,  U  that  as  the  absolute 
value  of  S  becomes  increasingly  smaller,  better  bounds  are  obtained  for  the 
net  radicat  flux.  The  values  of  the  radiant  flux  for  S  =  0,1  x  10"^  and 
«  =  -Ool  X  lG-3  aj^e  so  near  being  equal  to  each  other  that  the  curves  are 
coincident. 

It  vas  interesting  to  note  that  plots  of  the  maximum  and  minimum  re- 
siduals versus  S  are  shown  to  be  straight  lines. 

In  looking  for  methods  to  reduce  the  error  associated  with  the  N-*^ 
variation  in  the  assumed  form  of  the  power  series  solution,  it  is  recom- 
mended to  increase  the  value  of  N.  The  result  will  cause  the  bounds  to 

2U 
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approach  the  exact  solution  with  greater  accxxracy.  In  the  formolation  of  this 
problem;,  note  that  the  values  of  N  are  even  integers. 

The  surface  temperature  of  the  plate,  C7  (X),  vxas  asstaned  constant  for 

s 

the  purpose  of  con^aring  the  results  vxith  other  available  approximate  and 
exact  solutions.  With  the  present  program,  the  plate  temperature  may  be 
varied  at  different  positions  along  the  plate. 

CONCLUSIONS 

Although  complete  agreement  of  the  solutions  iJith  the  exact  solution 
has  not  been  obtained,  the  results  indicate  that  the  bounding  principle, 
applied  to  integral  equations,  gives  convergent  upper  and  lafer  bounds. 
The  procedure  can  be  extended  to  other  problems  which  are  described  by 
integral  equations.  Integration  of  equation  (l6)  is  a  relatively  simple 
operation  for  the  fu^'iction  Y  and  U,  However,  the  solution  becomes  more 
complex  as  the  geometry  of  the  problem  is  varied. 

It  is  believed  that  the  application  of  the  bounding  principles  will 
be  a  welcome  addition  to  the  field  of  radiation  space  technology  and  to 
the  solution  of  integral  equations  encountered  in  radiant  heat  transfer 
problems.  From  the  formulation  presented  in  this  paper,  it  is  now  possible 
to  deterr,iine  both  the  local  and  the  overall  heat  transfer  from  the  system. 

When  the  plates  which  exchange  heat  are  non-gray  radiators,  the  equations 
must  be  written  for  monochromatic  radiation. 

The  degree  of  accuracy  of  tlie  solution  obtained  by  boxmding  principles 
is  dependent  upon  the  efficiency  and  spead  of  present  day  high  speed  digital 
computing  facilities.  The  method  is  practical  and  affords  accurate  results. 
Moreover,  the  use  of  this  method  eliminates  the  need  for  "error  analysis". 
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APPENDIX  A 
Evaluation  of  the  Indefinite  Integral 
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EVALUATION  OF  TH2  INDEFINITE  INTEGRAL 


For  the  indefinite  integral  of  equation  (15),  Grobner  [fj  has  given 


the  integral  in  the  form 


where  N  ;^  2k  and   A  =    1    j    B  «=  -  (2N  ~  2k  ~  l)b  j 

CN  -  2v)a  (N  -  2k)a 


D  =  -  ^  -   l)g 
Tn  -  2k)a 


The  denominator  U  is  a  quadratic  equation  expressed  as 

Z/=  a.  /^v-  26/^  c  (26) 

For  the  radiant  flux  problem,  between  tv:o  piirallel  plates 


and  upon  ejqjanding 


(27) 


J^     ^^-  ZXY   ^    IX^I-H')  (27a) 


^ 


[.  30 

Thus  the  quantities  a,  b,  and  c  are  a  =  1,  b  =  -X,  and  c  =  X^  +  H^, 
The  value  of  k  in  equation  (25)  is  k  =  1  so  that  U^^'^s  r>  u*. 

RecrOl  that  the  requirement  N  7^  2k  must  be  satisfied  for  equation  (25). 
Since  k  =  1,  then  it  rToHoi-rs  that  K  r   2»  The  efiect  of  this  requirenient 
indicates  that  equation  (25)  can  be  used  for  any  value  of  N  in  the  pavter 
series  expansion  of  the  dur-jy  variable  Y  in  equation  (7a)  except  at  the 
value  N  =  2.  This  author  elected  to  use  the  integral  equation  for  values 
of  N  where  N  >  2  and  to  evaluate  the  integral  at  N  »  0,  1  and  2  by  a  standard 
table  of  integrals,  [j]   For  N  =  0,  the  integral  equation  becomes 


d^     =     /    ay^6      =   ^-^ 


ly^y^  ac-b^       j'/^  //^z/'/^ 


For  N  =  1 


/^/  =   -/   fZ-^c    ^   xY'X^-H 


U^/i  ac-  o"^       ijVz  y^  u'^^ 


Finally^  N  =  2  results  in 


/V/=  l\d'/     -t  /     (2b^-ac)Z  4.  be 


U^v^        aj  u^j^      a((xc-h^)  j'/2. 


(28) 


(29) 


(30) 
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d^  -     /     In  [cY-^b   -^  U 


lit 


(31) 


St;ibstituting  equation  (31)  into  equation  (30)  yields 


/V/  =  In  CZ-X-^^'^)  ^  xy-X'-XU'-  W^     (30a) 


APPENDIX  B 


Listing  of  Fortran  Program  to 
Solve  for  Coefficients,  A^. 


.  C  *. 
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KSU  14 IC  CC^.PUTING  CEKTER    .   ' 

BOP   CIEGC  RATH 
C      CCEFFICIENT  MATRIX  FCR  RADIATION  PROBLEM 
IC  CIPENSICN  B(16,I7) 

1  FCRNAT{5E16.8) 

2  FCRNAT (1H0,5EI6.8) 

3  FCR^AT(I3) 

I*;  REAC  INPUT  TAPE  5,3,« 

If-    REAC  INPUT  TAPE  5,1,CELX 

17  REAC  INPUT  TAPE  5,1,XG 

IP  REAC  INPUT  TAPE  5,1, XT 

1<=  REAC  INPUT  TAPE  5,1,H  ^ 

2C    REAC  INPUT  TAPE  5,1,RH0  '    .    ' 

2?  REAC  INPUT  TAPE  5  , 1 ,  ZAC  ,  ZCORR 

23  READ  INPUT  TAPE  5,1,T 

24  REAC  INPUT  TAPE  5,3,MN 

25  REAC  INPUT  TAPE  5,3,NC£LTA 
GIN=NIN 
DELX1=1.0/GIN 

NIN1  =  MN  +  1 
^2=(M/2}+l 
M3=N2+1 
LN  =  F2 

DC  60  MID=1,NDELTA 
21  REAC  INPUT  TAPE  5,1, DELTA 
9  FCRNAT(1H0,5X,2CH  DELTA  IS  NOW  SET  TO, £16. 8) 
WRITE  OUTPUT  TAPE  6, 9, DELTA 

5  FCRNAT(lH0,7X,AhC£LX,13X,2HXQ,15X,2HXT,14X,lHH,14X,3HRH0) 

WRITE  OUTPUT  TAPE  6,5 

WRITE  OUTPUT  TAPE  6 ,2 , DELX, XQ, XT,H ,RHC 

6  FCR^AT{1HO,7X,5HCELTA,11X,3H2AC,13X,5HZCORR,13X,1HT) 
WRITE  OUTPUT  TAPE  6,6 

WRITE  OUTPUT  TAPE  6 ,2 , DELTA, ZAC,ZCORR, T 

7  FCRNAT(1H0,7X,1HK,5X,3HMN,5X,6HNDELTA) 
WRITE  OUTPUT  TAPE  6,7 

8  FCRNAT(1H0,5X, I3,5X,I3,6X,I3) 
WRITE  OUTPUT  TAPE  6  ,8,  N,  ,MN,  NDELTA 
X  =  XC 

D£LX=2.»(XT-XQ)/CM 
DC  50  1=1, M2 
B{ I  ,N3)=DELTA+T««4 
N  =  2 

UPO=SGRTF( (l.-X)«(l.-X)4H«H) 
UNE=SCRTF( (l.+X)«{l.+X)+H«H) 
Z=(1.-X)/(H*H»UPC)+(1.+X)/(H«H*UNE) 
B(  I  ,l)=l.-(RHC«H«H/2.)«2 
Ll=l 

ZN2={X-X«X-H*H)/{h*H«UPC)+(X+X«X+H»H)/(H«H*UNE) 

ZM=LCGF{1.-X+UFC)+(X«X-X»X«X-X*H*H-H»H)/(H*H«UP0)-L0GF(-1.-X+UNE) 
l+(X«X+X«X«X+X»H»h-H«h)/(H«H*UNE) 
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B( I ,2)=X«X-(RHC*H«H/2.  )«2M 
L2  =  2 

CC    AO    J=3,M2 
DC    30    K=l,2 
N  =  N  +  1 
R  =  N 

ZN=l./( (R-2.)«UPC)+(2.«B-3.)*X»ZNl/(R-2. ) + ( ( R+ 1. ) / ( R-2 . ) )«{X»X+H«H 
l)»ZN2-( {-l.)»« (N-1) )/( (R-2.)«UNE)-( ( 2. »R-3 . ) / ( R-2. ) ) •X*2N1- ( (R+1.) 
2/(R-2.))«(X«X+H«H)*2N2 
ZN2  =  ZM 
3C    ZN1=ZN 

4C    B(  I  ,J)=X««N-(RHC«h»H/2. )*ZN 
5C    X=X+CELX 

MATRIX  SOLUTION  FCR  IBM  1410   2/64   APPL-RATHFON 
50?  FCRNAT(2H0  ,E16.8,2I3) 
5C3  DIMENSION  8BB ( 16 , 17 ) , A BA ( 16 ) , BSTAR ( 16 ) 

505  FCRMAT(2H0  ,2E16.8,I3) 

506  FCRMAT(2H0$,E16.8) 
514  FCRfAT(2H0  ,5£16,8) 
51?  ZK=9.CE10 

LM  =  LN  +  1 

DC  508  LI=1,LN 

ABA(LI )=0.0 

508  BSTAR(LI)=0.0 
ZTRV=0.0 
ZFA=0.0 

ZFAC=1.0     ■   ' 
ZZFAC=1.0 

DC  509  1=1, LN 
ZLL  =  I 

509  ZZFAC=ZZFAC»ZLL 
51C  ZZ=C.O 

CC  547  LI=1,LN  V    . 

LIR=LN+1-LI 

IF(B(LIR,LIR) )523,521,521 
521  ZY=E{LIR,LIR)-ZAC 
GC  TC  524 

523  ZY=-l.»e (LIR,LIR)-ZAC 

524  IF(ZY)525,525,547 
52'?  IF(LN-LI)526,526,527 

526  LI  1  =  2 
ZZ=1.0 

GC  TC  528  .      --  :    '^ 

527  LI1  =  LI  +  1  -    .     ■-• 
52P  CONTINUE 

DC  537  LK=LI1,LN  ,-■ 

IF(ZZ)529,529,530        '         -  [  ■; 
529  LK1=LN+1-LK 

GC  TO  531 
53C  LK1=LK 
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531  IF(6(LK1,LIR))532,533,533 

532  2Y=-1.*B(LK1,LIR)-ZAC 
GC  TC  534 

533  ZY=E(LK1,LIR)-2AC 

534  IF{ZY)535,535,541 
53*5  IF(  ZK-ZY)537,537,536 
53f  ZK=ZY 

537  CCNTINUE 

IF(ZZ)538, 538,540 
53P  LI1=LIR+1 

ZZ  =  1. 

P  P   T  P  S  ^  fl 

53<5  FCRNAT(1H0,23H  ALL  ELEMENTS  IN  COLUMN,  13 ,21H  ARE  SMALLER  THAN  ZAC) 

540  ZK=ZK+ZAC 

WRITE  OUTPUT  TAPE  6,539,LIR 
WRITE  OUTPUT  TAPE  6, 514, ZK, ZAC 
ZAC=ZAC/10. 
GC  TC  510 
C      REACY  TC  INTERCHANGE  TWC  ROWS 

541  DC  542  LJ=1,LN1 
ZTE=B(LIR,LJ) 
B(LIR,LJ)=B(LK1,LJ) 

542  B(LK1,LJ)=ZTE 
ZFA=ZFA+1.0 
ZFAC=ZFAC»2FA 

IF (2ZFAC-ZFAC) 544,544,545 

543  FCRMAT(1H0,17H  ZAC  CECREASED  T0,E16.8,25H  IN  LIEU  OF  REARRANGEMENT 

1) 

544  ZAC=ZAC/10. 

WRITE  OUTPUT  TAPE  6, 543, ZAC 

545  IF(ZZ)547,547,546  , 
54  6  GC  TC  510 

547  CCNTINUE 

IF(SENSE  SWITCH  1)548,550  .   " 

548  DC  549  LI=1,LN 
DC  549  LJ=1,LN1 

54<5  WRITE  OUTPUT  TAPE  6  ,502  ,B  (LI  ,  LJ) ,  LI  ,LJ 

55C  CCNTINUE  .     ' 

DC  551  LI=1,LN 

DC  551  LJ=1,LN1 

551  BBB(LI,LJ)=B(LI,LJ) 

C      MACHINE  IS  DONE  REARRANGING  AND  READY  TO  SOLVE  THE  MATRIX 
511  DC  555  LI=1,LN 
ZELC  =  B(LI,Ln 
DC  552  LJ=LI,LN1 

552  B(LI,LJ)=B(LI,LJ)/ZELD 
DC  555  LK=1,LN 
IF(LK-LI)553,555,553 

553  ZELS=8(LK,LI) 
DC  554  LJ=1,LNI 
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55A  B{LK,LJ)=B(LK,LJ)-Z£LS»B(LI,LJ) 
55«  CCMINUE 

ZTRV=ZTRY+1.0 
C      MATRIX  IS  SCLVEC 

IF(SENSE  SWITCH  2)556,558 

556  WRITE  OUTPUT  TAPE  6,51A,2TRY 
DC  557  LI=1,LN 

557  WRITE  OUTPUT  TAPE  6 ,502  ,E (LI , LNl) , LI 
55P  DC  559  LI=1,LN 

55<5  ABA(LI)  =  B(LI,LNl)+ABA(LI) 

CONTINUE 

CC  560  LI=1,LN 

DC  560  LJ=1,LN 

BSTAR{LI )=ABA(LJ)*BEB(LI,LJ)+BSTAR(LI ) 
56C  8(LI,LJ)=BBB(LI,LJ) 

CC  561  LI=1,LN 

B(LI,LNl)=0BB(LI,LNl)-BSTAR(Ln 

561  BSTAR(LI)=0.0 

IF (ZTRY-ZCORR) 562,565,565 

562  ZZE=C.0 

CC  563  LI=1,LN 

563  ZZE=ZZE+B(LI,LN1) 
IF(ZZ£)564,565,564 

56A  GC  TO  511 

565  WRITE  OUTPUT  TAPE  6  ,514  ,  2TRY , ZCORR 
CC  566  LI=l,LN 

566  WRITE  OUTPUT  TAPE  6 ,505 , ABA ( LI ) , B { L I , LNl ) , LI 
X  =  0.0 

RNAX=-1.0E+06 

RNIN=l.CE+06 

DC    150    L=1,NIN1 

N  =  2 

UPC=SCRTF( {l.-X)«(l.-X)4H«H) 

UNE=SCRTF( {l.+X)«(l.+X)+H«H) 

Z={1.-X)/IH«H«UPC)+{1.  +  X)/(H«H«UNE)  '■ 

RES=ABA(1)*(1.-(RHC«H«H/2.)«Z) 

ZN2={X-X«X-H»H)/(H»H«UPC)+(X+X»X+H»H)/(H«H«UNE) 

ZM=LCGF{1.-X+UPC)+(X«X-X*X»X-X*H»H-H«H)/(H*H«UP0)-L06F(-1.-X+UNE) 
1+(X«X+X*X«X+X«H*H-H«H)/(H«H*UNE) 

RES=RES+ABA(2)*(X«X-(RHC«H«H/2.)»ZN1) 

CC    140    J=3,N2 

DC    130    K=l,2 

N  =  N-H  ' 

R  =  N 

ZN=l./( (R-2.)»UPC)+{2.*R-3.)*X«ZNl/(R-2.)+((R+l.)/(R-2.))«(X»X+H»H 
l)«ZN2-({-l.)««(N-l) )/( (R-2.)*UNE)-( ( 2.*R-3. ) / {R-2. ) ) *X»ZN1- ( (R+1. ) 
2/{R-2.))«(X«X+H*H)«ZN2 

ZN2=ZN1 
13C    ZM  =  ZN 
140    RES=RES+ABA{J)«(X««N-(RHC*H*H/2.)«ZN) 
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RES=RES-T«*4 

IF  (RES-RMAX) 102,102,101 
ICl    RNAX=RES 

XNAX=X 
1C2    IF{RNIN-RES)10A,1CA,103 
1C3    RMN  =  RES 

XNIN=X 

IFtSENSE  SWITCH  3)10^,150 
IC^  URITE  OUTPUT  TAPE  6, 514, RES, X 
15C  X=X+CELX1 
1C5  FCR^AT(1HC,6X,4HR^AX,12X,4HXMAX,12X,AHRMIN,12X,4HXMIN) 

WRITE  OUTPUT  TAPE  6,105 

WRITE  OUTPUT  TAPE  6  ,2 , RN AX ,XMAX,RMIN,XMIN 

IF(SENSE  SWITCH.  A)1C6, 60 
1C6  DC  107  LI=1,LN 
107  WRITE  OUTPUT  TAPE  6  ,506 ,AeA( LI ) 

60  CONTINUE 
lOS  CALL  EXIT 

STCP 

END      ■ 
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Listing  of  Fortran  Program  to 
Solve  for  Bounds 
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BOP   ECUNO  RATH 
C      DETERNINATICN  CF  UPPER  ANC  LOWER  BOUNDS  OF  BETA(X) 
DIN.EKSICN  AL0W(16),AUP(16) 

1  FCRKATCSEld.a) 

2  FCRNAT(1H0,5E16.8) 

3  FCRf'ATda)  . 

15  READ  INPUT  TAPE  5,3,M 

16  BEAC  INPUT  TAPE  5,1,CELX 

17  READ  INPUT  TAPE  5,1,XQ 
IP  REAC  INPUT  TAPE  5,1, XT 
1<;  REAC  INPUT  TAPE  5,1,H 
2C  REAC  INPUT  TAPE  5,1,RH0 

2?  REAC  INPUT  TAPE  5  , 1  ,ZAC  ,  ZCORR  ' 
23  REAC  INPUT  TAPE  5,1,T 
2A  REAC  INPUT  TAPE  5,3,NIN 
25  REAC  INPUT  TAPE  5,3,NDELTA 

N2={M/2)+l 

DC  9  NDELTA=1,2 
21  REAC  INPUT  TAPE  5,1, DELTA 

5  FCRKAT(1H0,7X,AHCELX,13X,2HXQ,15X,2HXT,IAX,1HH,1AX,3HRH0) 

WRITE  OUTPUT  TAPE  6,5 

WRITE  OUTPUT  TAPE  6  ,2 , DELX , XQ, XT,H,RHO 

6  FCRNAT(1H0,7X,5HCELTA,11X,3HZAC,13X,5HZC0RR,13X,1HT) 
WRITE  OUTPUT  TAPE  6,6 

WRITE  OUTPUT  TAPE  6  ,2 , DELTA, ZAC,ZCORR,T 

7  FCR>'AT(1H0,7X,1HN,5X,3HMN,5X,6HNDELTA) 
WRITE  OUTPUT  TAPE  6,7 

8  FCR^AT{1H0,5X, I3,5X,I3,6X,I3) 

«3  WRITE  OUTPUT  TAPE  6  ,8,  M  ,NIN,  NDELTA 

BEAC  INPUT  TAPE  5,1,DARX 

REAC  INPUT  TAPE  5, 3, MEL 

DC  30  N=1,H2 

REAC  INPUT  TAPE  5,1,AL0W(N) 
3C  WRITE  0UTPUTTAPE6,2,AL0W{N) 

DC  AO  N=1,M2 

REAC  INPUT  TAPE  5,1,AUP{N) 
4C  WRITE  OUTPUT  TAPE  6,2,AUP(N) 

IC  FCRNAT(1H0,2H  X  ,  5X , IIHLCWER  BOUND, 5X, IIHUPPER  BOUND, 5X, lOHMEAN  VAL 
1UE,5X,17HNAX  PERCENT  ERROR) 

WRITE  OUTPUT  TAPE  6,10 

X  =  0.0 

DC  50  1=1, MEL 

BETAL=0.0 

BETAU=0.0 

BETAL=AL0W{1)+BETAL 

BETAU=AUP(1)+BETAU 

DC  60  N=2,M2 

BETAL=AL0W(N)«X*«{2*N-2)+BETAL 
6C  BETAU=AUP(N)«X#«(2«N-2)+BETAU 

B£TA«=(BETAU+BETAL)/2.0 
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TAXER»( (BETAU-BETAL)/(2.0«BETAL))»100.0 
11  FCR^AT(IHC,F4.2,^E16.8) 

WRITE  OUTPUT  TAPE  6 , 11 t X,BETAL, BETAU, BETAM.TAXER 

5C  X»X+DARX 
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A  method  of  bounding  principles  is  presented  for  solving  the  integral 
equation  which  arises  in  determining  the  combined  radiation  flux  between 
two  non-black,  finite,  parallel  plates.  Convergent  upper  and  lower  bounds 
are  obtained  in  order  to  approximate  the  solution  to  the  problem.  The 
method  emphasizes  the  attainment  of  an  approximate  solution  of  known  accuracy, 
is  systematic,  and  is  well  adapted  to  analysis  on  high  speed  digital  computers, 

Agreement  with  existing  solutions  has  not  been  achieved  and  hence,  addi- 
tional study  is  needed,     ■ 
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